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Abstract 

This paper is concerned with coherent quantum control design for translation invariant networks of identical quantum 
stochastic systems subjected to external quantum noise. The network is modelled as an open quantum harmonic oscillator and 
is governed by a set of linear quantum stochastic differential equations. The dynamic variables of this quantum plant satisfy 
the canonical commutation relations. Similar large-scale systems can be found, for example, in quantum metamaterials and 
optical lattices. The problem under consideration is to design a stabilizing decentralized coherent quantum controller in the 
form of another translation invariant quantum system, directly coupled to the plant, so as to minimize a weighted mean square 
functional of the dynamic variables of the interconnected networks. We consider this problem in the thermodynamic limit of 
infinite network size and present first-order necessary conditions for optimality of the controller. 

I. INTRODUCTION 

Currently emerging technologies open up opportunities to synthesize artificial optical media known as quantum 
metamaterials; see, for example, [17], [12], [29], [16], [27]. These large-scale quantum networks with are engineered from 
complex unit cells and are effectively homogeneous (in the sense of translational symmetries) on the scale of relevant 
wavelengths, for example, in the microwave range. In natural solids, the quantum energy level configurations of the constituent 
atoms or molecules specifies the optical behaviour of the material. 

In contrast, the controllable resonant characteristics of the building elements of quantum metamaterials, such as the 
Josephson devices or optical cavities [17], [12], [16], determine the electromagnetic response of the quantum metamaterials. 
The quantum metamaterials are considered to be a promising approach to the implementation of quantum computer elements 
which can maintain quantum coherence over many cycles of their internal evolution [16], [26]. 

A framework for modelling and analysis of a wide range of open quantum systems, including those arising in quantum 
metamaterials, is provided by quantum stochastic differential equations (QSDEs) [8], [13]. In QSDEs, the environment is 
modelled as a heat bath of external fields acting on a boson Eock space [13]. In particular, linear QSDEs represent the 
Heisenberg evolution of pairs of conjugate operators in a multi-mode open quantum harmonic oscillator which is coupled 
to the external bosonic fields. This framework also allows for robust stability analysis for certain classes of perturbed open 
quantum systems which have been addressed, for example, in [15], [19]. 

The analysis of large-scale quantum networks can be effectively reduced in the case when they are organized as a 
translation invariant interconnection of identical elements with periodic boundary conditions (PBCs); see, for example, [16], 
[24], [22]. The PBCs rely on negligibility of boundary effects in such a network consisting of a sufficiently large number 
of subsystems. This technique is used for lattice models of interacting particle systems in statistical physics (for example, 
in the Ising model of ferromagnetism [11]). 

Coherent quantum feedback control is aimed at achieving robust stability and robust performance through measurement- 
free interconnection of quantum systems [9], [20]. In this approach, the controller is another quantum system which is 
coupled to the quantum plant, for example, through a bilateral energy interconnection, known as direct coupling [28]. In 
comparison with the more traditional measurement-based feedback control techniques, coherent control benefits from the 
preservation of quantum coherence within the network. 

In this paper, both the quantum plant and the coherent quantum controller are modelled as large fragments of translation 
invariant networks of linear quantum stochastic systems endowed with the PBCs. The nodes of the networks are directly 
coupled to each other within a finite interaction range. This interconnection is governed by linear QSDEs based on the 
Hamiltonian and coupling parametrization of the corresponding multi-mode open quantum harmonic oscillator whose 
dynamic variables satisfy the canonical commutation relations (CCRs). Eollowing a similar approach used in the classical 
control theory, we employ spatial Eourier transforms in order to obtain a more tractable representation of the dynamics of 
the quantum feedback network; see, for example, [25], [2], [24], [22] and the references therein. 

We consider a weighted mean square performance index for the stable plant-controller network in the thermodynamic 
limit (when the network size goes to infinity). In this framework, the decentralized control design problem is formulated as 
the minimization of the cost functional over the parameters of the controller and its coupling with the plant. By calculating 
the Erechet derivatives of the cost functional, we obtain first-order necessary conditions for optimality of the controller. In 
comparison with previous results on coherent quantum LQG control [23] using Gramians and related algebraic Lyapunov 
equations (ALEs), the optimality conditions developed in the present paper employ spatial spectral densities, which reflects 
the translation invariant nature of the underlying problem. 

This work is supported by the Australian Research Council. The authors are with UNSW Canben'a, ACT 2600, Australia. arash_kho@hotmail. com, 
igor.g.vladimirov@gmail.com, i.r.petersen@gmail.com. 


II. NOTATION 


Unless specified otherwise, vectors are organized as columns, and the transpose (•)^ acts on matrices with operator-valued 
entries as if the latter were scalars. For a vector X of operators Xi,...,Xr and a vector Y of operators Ti,..., 7^, the commutator 
matrix is defined as an (r x s)-matrix :=XY^ — whose ( 7 ,A;)th entry is the commutator —Tj-Xy 

of the operators Xj and T^. Also, (•)^ := ((•)**)^ denotes the transpose of the entry-wise operator adjoint {■)*. In application 
to complex matrices, reduces to the complex conjugate transpose (•)* := ((•))^- Furthermore, Sr and Ar denote the 
subspaces of real symmetric and real antisymmetric matrices of order r, respectively and A 2 is spanned by the matrix 
J := g . Also, i := denotes the imaginary unit, Ir is the identity matrix of order r, positive (semi-) definiteness 
of matrices is denoted by {'^) and 0 is the tensor product of spaces or operators (for example, the Kronecker product 
of matrices). The spectral radius of a matrix M is denoted by r(M). The adjoints and self-adjointness of linear operators 
acting on matrices is understood in the sense of the Frobenius inner product {M,N) := Tr{M*N) of real or complex matrices. 
The Kronecker delta is denoted by 5jk, and U := {z S C : |z| = 1} is the unit circle in the complex plain. The complex 
residue of a function / about a point zo is denoted by Rss,z=zof(.^)- Also, := Tr(pi^) denotes the quantum expectation 
of a quantum variable ^ (or a matrix of such variables) over a density operator p which specifies the underlying quantum 
state. For matrices of quantum variables, the expectation is evaluated entry-wise. 


III. LINEAR QUANTUM STOCHASTIC SYSTEMS 

We consider a quantum stochastic system interacting with external boson fields [8], [13]. The system has N subsystems 
with associated 2n-dimensional vectors Xq, ... ,Xn-i of dynamic variables which satisfy the CCRs 

r Yo 

[X,X'^]=2i&, 0:=/a,0 0, Y:= : 

Yw-i. 


( 1 ) 


Here, 0 is a block diagonal joint CCR matrix, where © G A 2 „ is a nonsingular matrix. The system variables evolve in time 
according to the QSDE r wh -1 

( 2 ) 


dX=(i[H,X]-^BJB'^0^^xyt+BdW, W : = 


Wn- 


Here, Wb,...,lUv_i are 2m-dimensional vectors of quantum Wiener processes with a positive semi-definite Ito matrix Q. G 

El2m: 

dWjdW^ = 5jkQ.dt , a:=l2„ + iJ, 7: = 4 0 J. (3) 

Accordingly, J := /Ar07 is a block diagonal matrix, and the matrix B G ]R2«wx2mV jjj (|2]) is related to a matrix M G ]R 2 mA'x 2 «w 
of linear dependence of the system-field coupling operators on the system variables by 

B:=20M^. (4) 


The term —jBJB^O ^X in the drift of the QSDE (|2| is associated with the system-field interaction and results from 
evaluating the Gorini-Kossakowski-Sudarshan-Lindblad decoherence superoperator [6], [10] at the system variables. Also, 
H is the Hamiltonian which describes the self-energy of the system and is usually represented as a function of the system 
variables. In the case of an open quantum harmonic oscillator [4], [5], the Hamiltonian H is a quadratic function of the 
system variables 

1 1 

H := -X^RX = - ^ xjRjkXk, (5) 

^ ^ j.k=Q 


where the energy matrix R := {Rjk)o^j,k<N S S2nN is formed from blocks Rj^ = RJj G By substituting Q into (|2]) 

and using the CCRs ([T]l, it follows that the QSDE takes the form of a linear QSDE 

dX = AXdt -f BdW, (6) 

where the system matrix A G ]R 2 nA'x 2 «w given by 

A:=20R--BJB^0 \ (7) 

2 

We consider a 2^A-dimensional vector Y of q'-dimensional output fields Yq,. .. ,Yt^_\ which is evolved by the linear QSDE 


dT = CXdt A- DdW, 


■ To ■ 

Jn-i. 


( 8 ) 
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with C £ '^}^Nx 2 nN 2) g ^ satisfies the non-demolition condition [3] with respect to the dynamic variables 

in the sense that 

[X(f),F(5f]=0, t^s. (9) 

Due to this non-demolition nature, the output fields can be interpreted as ideal measurements over the open quantum system. 
The condition (|9]) implies an algebraic relation between C and D in ®: 

0C^^BJD^ (10) 

Furthermore, the matrix DJD^ is an appropriate submatrix of J from Q (so that q ^ m). 

IV. LINEAR QUANTUM STOCHASTIC NETWORK WITH PERIODIC BOUNDARY CONDITIONS 

Suppose the open quantum harmonic oscillator of the previous section represents a fragment of a translation invariant 
network which is organised as a one-dimensional chain of identical linear quantum stochastic systems numbered hyk — 
0,... ,A — 1. Each node of the network interacts with the corresponding external boson field, and hence, the joint network-field 
coupling matrix M in (|4]l and the feedthrough matrix E) in ® are block diagonal; 

M-=In®M, D=In®D, (11) 

where M G R 2 mx 2 n g Hence, the matrix B is block diagonal and so also is the matrix C in view of (fTOl i. 

(fTTIl and nonsingularity of the CCR matrix 0 in ([T]i: 

B=In®B, C = In®C, B:=20M'^, C:=2DJM. (12) 

The nodes in the network are directly coupled to each other within a finite interaction range d. The fragment of the chain 
is assumed to be large enough in the sense that N > 2d, and is endowed with the PBCs, thus having a ring topology. A 
particular case of nearest neighbour interaction (when d =\) is depicted in Pig. [T] In the case of an arbitrary interaction 



Yk+l i - F i - > F -^ 



Fig. 1. A finite fragment of the translation invariant network of open quantum systems with direct coupling between the nearest neighbours. Also shown 
are the input and output fields of the nodes of the network. 


range d^ I, the Hamiltonian // in (|5]l is completely specified by matrices Ri = ^ with £ = 0,±1,.. .,±d (so 

that Rq £ S 2 n) as 

^N-i d 

^ = T L ^7 L • (13) 

^ ;=0 ^ t=-d ' 

Here, j — i is computed modulo N in accordance with the PBCsHl and hence, the corresponding matrix R is block circulant. 
The matrix Rq € S 2 „ specifies the free Hamiltonian for each node, while R±s describe the energy coupling between the nodes 
which are at a distance s= \ . ,d from each other (with the more distant nodes in the network not being directly coupled). 

In view of the block diagonal structure of the matrices 0 and M in ([T]i and (fTTT i. it follows from (|4|i, (|7]i and (fTJt that the 
QSDE (|6]l is representable as a set of coupled QSDEs for the dynamic variables of the nodes of the network; 

d 

AXj = A^V„„d(2„£,iV)df +BdW;, ; = 0,... ,A- 1, (14) 

e=-d 


where use is made of (fT^ . and the matrices A_j,... ,Aq £ are given by 


r 2@Rq-\BJB^@-^ iff = 0 

I 2 @Rt if 0 <\e\^d 


( 15 ) 


^with the MATLAB function being used instead of the standard modular arithmetic notation in order to avoid confusion in the subscripts 
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Furthermore, the output fields Yq,... jFat-i in (l8]l, associated with the nodes of the network, evolve according to the QSDEs 


AYj = CXjdt + Dmj, j = Q,...,N-l. (16) 

where DJD^ is a submatrix of J in view of the block diagonal structure of J and D. Therefore, the network, governed by 
(fl4li and (fTfil l. has a block circulant structure. 

V. INTERCONNECTED PLANT AND CONTROLLER NETWORKS 

Consider an interconnection of two translation invariant networks, similar to the one described in Section IIVI where 
one of the networks is a quantum plant while the other is interpreted as a decentralized coherent quantum controller. We 
consider the case when the plant and controller networks are of equal size N. Eor example, Eig. |2] provides a layout of the 
plant-controller network with nearest neighbour direct coupling. In what follows, the operators and parameters of the plant 


Y\,k-\ 



Yi,k+\ 


Fig. 2. Finite fragments of two translation invariant quantum networks of equal size with nearest neighbour direct coupling. The 1st network is a plant 
and the 2nd is a controller. Also shown are the input and output fields of the nodes of the networks. 


and controller networks will be indicated by subscripts 1 and 2, respectively. In particular, X\ j and Xjj denote the vectors 
of dynamic variables for the y'th nodes of the plant and the controller, respectively, with the corresponding dimensions 2ni 
and 2n2- The total Hamiltonian of the plant-controller interconnection is 


H:^Hi+H2 + Hn. 


Here, in view of (fOT l. the plant Hamiltonian Hi and the controller Hamiltonian Hi are given by 

^ i^kj ^ ^k/^k,mod{j—i,N)p ^ 

^ j=o e=-dt 

where S R 2 ntx 2 nj.^ £ = o,±l,.. .,±dk. Also, Hu is the interaction Hamiltonian; 


7V-1 , £/ _ 

Hn := E fe E 

J=o 




i=-d 


(17) 


(18) 


(19) 


where d denotes the range of direct coupling between the plant and controller nodes (with N > 2max{di,d2,d)), and 
R( G ]R 2 nix 2«2 ^j-g j-jjg plant-controller coupling matrices, with |f| ^ d. It is assumed that the plant variables commute with 
the controller variables: [Vi = 0 for all 0 ^ j,k < N. Hence, Hu in ([19]) is indeed a self-adjoint operator. Due to the 

structure of the total Hamiltonian, the plant-controller network is a block circulant system [25]. By substituting H specified 
by (I17l) - (ll91 l into the QSDE (|2]i, it follows that the plant and controller variables are governed by a set of coupled QSDEs: 





d^2.; 


+ BidWij, 

l=-d e=-d2 


+ B2dW2j 


( 20 ) 


( 21 ) 
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for all 0 ^ j < N, where, similarly to (fTSl i and (fTSl) . the matrices of coefficients are given by 

Au-=^ 20 ^^i,o - if f = 0 


2&kRk,i if 0 < \i\ < dk ’ 

Aij := 2&iRe, A 2 ,f := 202 ,^, \i\ < d, 


Bk:=2QkMl. 


( 22 ) 

(23) 

(24) 


The coupled QSDEs (l20l) - (l2ni . which are associated with nodes of the plant and the controller, can be studied in the spatial 
frequency domain. Similarly to [24], we will use the discrete Fourier transforms (DFTs) of the quantum processes Xkj and 
Wk.j over the spatial subscripts j = Q,...,N — 1 for k— 1,2: 


■= nAt) ■= 

j=0 j=0 

which are considered for z G Uat from the set of Mh roots of unity 

Uw:={eT : j = 0,... ,N - l}. 

Then the corresponding augmented quantum process is evolved in time by the QSDF 


^2,z 


n,z 

n,z 


djr, = jr.df + mw,, := 

where the matrices Aj and are defined in terms of the network parameters in (l22T i- (l24l i as 




, ^ := 

Bi o' 




0 B 2 


(25) 


(26) 


(27) 


(28) 


see [22] for more details. Here, use is made of (I25b together with the PBCs and the well-known properties of DFTs due to 
the assumption that z G Uat. As discussed in [24], it follows from the CCRs ([T]i that 


j,k=0 

iV-1 


k=0 


= 2i@YQ' = Xj:= Ij 


(29) 


for all z,v G Utv, where 0 := q' q 
a similar reasoning, Q implies that 


is the common CCR matrix of the augmented vectors Xj for all 7 = 0,... ,A — 1. By 




a- 


0 


0 

0.2 ’ 


(30) 


where Oi and Oj are the quantum ltd matrices of the plant and controller Wiener processes VTi ,■ and VT 2 . 7 . Assuming the 
network size N to be fixed, we denote the matrix of second-order cross-moments of the quantum processes and in 
(ED by 

(31) 


Here and in what follows, the quantum expectation E( ) is over the tensor product p (U of the initial state U5 of the 
plant-controller network and the vacuum state u of the external fields. 

Lemma 1: Suppose the matrix in (i28l l is Hurwitz for any z G Uat from (i26l l. Then for all z,v G Uat, the matrix in (i3Tl l 
has the limit .^z.A°°) •“ linir^+<=o'^.v(0 =Nd^v^z’ where the matrix 0 is a unique solution of the AFF 

+ S^zK + = 0. (32) 

The proof of Femma [T] can be found in [22] and is based on linearity of the QSDF ( iZTi i. the relation ( [30b and the fact 
that the forward increments of the quantum Wiener process in the vacuum state are uncorrelated with the adapted processes. 
The matrix-valued function U 9 z 1 —in (l32l i (which is well-defined on the unit circle under a stronger stability condition 
formulated in the next section) is the spatial spectral density [24] of the closed-loop system. This density encodes the 
covariance structure of the plant-controller variables in the invariant Gaussian quantum state [14] in the limit of infinite 
network size N ^ -l-oo. In this sense, ,5^ is a network counterpart to the controllability Gramian [1] of the system. 
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VI. OPTIMAL COHERENT QUANTUM CONTROL PROBLEM 

In view of the block circulant structure of the plant-controller network of Section |Vl this closed-loop system is stable for 
any network size V ^ 1 if and only if the matrix in (|28]) is Hurwitz for all z S U, that is. 


maxr(e-^Q < 1. 

ZGU 


(33) 


Here, use is also made of continuity of with respect to z G U and the property that U 7 v=i ^ dense subset of U. In 

what follows, the plant-controller network is called stable if it satisfies (l3^ . In this case, the decentralized coherent quantum 
controller is referred to as a stabilizing controller. As a performance criterion for such controllers, we will use a steady-state 
weighted mean square cost functional which has to be minimized in the limit of infinite network size. Eor simplicity, it is 
assumed that each node of the controller is only coupled to one node of the plant, that is, d = 0. The performance cost 
functional is given by 




^ lim E Y, Xjitfoj^kXkit), 


(34) 


where Xj are auxiliary quantum processes associated with Xj from ( |29] | by 


Xj{t):^EXj, 


hni 0 

0 Ro ’ 


(35) 


for y = 0,...,A — 1, and Rq is the coupling matrix which, in accordance with ( fT9l l, specifies the interaction Hamiltonian 
Hi 2 = in the case d = 0. Moreover, Ok is a given -valued sequence which satisfies a^k = O’/ for all 

integers k and specifies a real symmetric block Toeplitz weighting matrix {<yj-k)o^j,k<N- The block Toeplitz structure of the 
weighting matrix in 0^ corresponds to the translation invariance of the quantum network being considered. A matrix-valued 
map U 9 z = E*, defined as the Eourier transform 


E, := Y (36) 

k= — oo 

describes the spectral density of the weighting sequence. In order to ensure the absolute convergence of the series in (l3^ . it 
is assumed that L/T-oc II 0*11 < which also makes E^ a continuous function of z. The fulfillment of the condition E^ )= 0 
for all z G U is necessary and sufficient for {(yj-k)o^j,k<N ^ 0 to hold for all A ^ 1; see, for example [7]. In this case, the 
sum on the right-hand side of (l34l) is a positive semi-definite operator, and hence, ^0. The cost functional (on in (IT^ 

resembles the classical LQG control performance index [1] which imposes a cost on the actuation signal. A similar approach 

can be used in order to “penalize” other variables of interest by an appropriate choice of the weighting matrices Gk in ( lT4l i 
and the matrix E in (iTSl l. Eor any stabilizing controller and any given (sufficiently large) network size N, Lemma [T] allows 
the cost functional in (l34l i can be computed as [22]: 

<^^ = 4 L (37) 

ZGUjv 

where the matrix is the unique solution of the ALE (l32l l. and T,n{z) '■= L|/i|< 7 v(l ~ Due to the uniform 

convergence limAr^+ooEAr(z) = E, to the spectral density (I36I) over z G U, the representation dJTl i leads to the following 
infinite network size limit of the mean square cost functional Sn in (l34l i: 

^ := lim 4/ = — / (E,, E5^,E'^) — 

2m Jv z 

(4-., = Rj (38) 

for any stabilizing controller in the sense of ( l3Tt : see [24], [22] for more details. The resulting cost functional (S' in (IMt 
corresponds to the thermodynamic limit of equilibrium statistical mechanics [18]. We will now consider a decentralized 
coherent quantum control problem which is formulated as the minimization 

S’ —>■ min (39) 

of the infinite network cost functional in (I38I) over the energy and coupling matrices of stabilizing controllers (with all the 
dimensions of individual nodes being fixed). A particular case of such a network with nearest neighbour interaction (when 
di = d 2 = I and = 0) is depicted in Eig. [3 With the matrix B 2 G ]R 2 n 2 x 2 m 2 jjj being fixed, the cost S in is a 
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<- > Fi < -^ Fi < - >■ Fi <- 


.\ 


Y 2 ,k \ 


i: 


\ Wi,t 

\ 




F 2 < - > F 2 <■ 

\W2,t-l 


-> ^2 <- 


\ 


\W 2 ,, 


Fig. 3. The infinite plant-controller network with nearest neighbour direct coupling. Also shown are the input and output fields. 


function of 


8 ■= (^2,0,^2,l) ■ ■ ■ ^R242,Rq) 

e §2^2 X (K2"2><2«2)</2 X K2«1X2«2 =. G (40) 

which parameterizes the controller matrices in (I22l l- (l2?t . The minimization of the cost functional S' in (|39] | is carried out 
over the set 

Go := S G : 22 ^^ in (1281) satisfies (1^ 1 (41) 

of those g which specify stabilizing quantum controllers (ISTT l for the quantum plant (|20] |. The set G on the right-hand side of 
(l40l i is endowed with the structure of a Hilbert space with the direct sum inner product {g,g') := {Ro,Rq) +L^io(^2/j^2 e)- 
In what follows, we will use an auxiliary spectral density dehned on the unit circle z G U as the unique solution of the 
ALE 


.e// = 0 (42) 

for a stabilizing controller, where use is made of (ITST i and (l36l l. The function corresponds to the observability Gramian 
[1]. We will also use a network counterpart to the Hankelian [23]: 


^ := 


(43) 


where is given by ( l32l i. These and related matrices are partitioned into four blocks (•)j^ indexed by 1 < j,k 2 
according to their association with the plant and controller variables, and the Fourier parameter z will sometimes be omitted. 
The following theorem, which is formulated for the case d = 0, can be extended to arbitrary range of the plant-controller 
coupling. 

Theorem 1: Suppose the plant-controller coupling range in the network is = 0. Then necessary conditions of optimality 
for a stabilizing controller in the problem (|39] | are as follows: 


p^^^ Re(z^(@2.^22-^22Q2)) Q 

z=0 Z 


' = 0,. . . ,^f2, 


(44) 

(45) 


p^^^ Re(,^/i02-@iJf]2 + ;(^.£^z)22) Q 

z=0 Z 

Proof: The matrices R 24 = (l22l l influence the cost functional S in ( l38b only through the matrix in 

which depends on those matrices through 22 ^ in (1281) . Hence, the first variation of the integrand in (l38l) with respect to R 24 

= (46) 


for any z G U and £ = 0,... ,d 2 , provided the controller parameters in (l40l l satisfy g G Go in (HTt . Here, use is made of the 
ALEs (I 32 I 1 . (l42l) and the Hankelian from (l43T l. and 




0 0 
0 2 


J 02 5^2.0 iff = 0 

\@2{z-‘8R2,iF^dRl f if f > 0 • 


Substitution of (l47l l into (l46l l yields the following Frechet derivatives 




2Re( ©2 - © 2 ^ 2 ) if f = 0 
4Re(z^(,^*2©2 - © 2 .^ 2 )) if f > 0 ’ 


(47) 


(48) 
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where the symmetry of /?2,o is taken into account. Upon integration according to 

= 


, it follows from (l48T l that 


2Res,=0 


if f = 0 

4Res iff>0 


( 49 ) 


Now, the plant-controller coupling matrix Rq influences the cost functional S’ not only through the matrices (|2^ . (l28T l and 
the ALE (I 32 I 1 . but also through the matrix E in (l35T l. Therefore, by appropriately modifying (l46l l. it follows that 


Ey,E'^)=2Rc{{Mi,, + {i:,Ey„ 5^E)), 

which holds for any stabilizing controller (with g G Go), where 


5 b £/y = 2 


0 

©25^ 


0i5/?o 

0 


«o 


0 0 
0 SRq 


(50) 


(51) 


Substitution of dSB into (I 50 I 1 leads to the corresponding Frechet derivative of the integrand in (ISFt : ES^^E^) = 

4 Re{jS2*i@2-@i^i2 + j{^zE^z) 22 ), and hence, 

Re(J^*i©2 - ©1^2 + 3 {^zES^z) 22 ) 


S = 4Res ■ 
^0 z=0 


(52) 


The necessary conditions of optimality (l44l) and (|45]) can now be obtained by equating to zero the Frechet derivatives in 
(|4^ and (l52T i. respectively. ■ 

The optimality conditions (l44l i and (l45T l provide a set of algebraic equations for finding an optimal controller in the problem 
(I 39 I 1 . Although their solution is not yet available, the Frechet derivatives of the cost functional in (l49l l and (l5^ can be used, 
for example, in a gradient descent numerical algorithm for finding locally optimal stabilizing controllers, similar to [21]. 


VII. CONCLUSION 

We have considered a decentralized coherent quantum control problem for a quantum plant modelled as a large-scale 
one-dimensional translation invariant network governed by linear QSDEs with PBCs. The problem is to design a similar 
coherent quantum controller network, directly coupled to the plant, so as to minimize a weighted mean square functional in 
the thermodynamic limit. We have used DFTs in order to achieve a more tractable form of the problem in the spatial frequency 
domain. This reformulation has allowed necessary conditions of optimality to be obtained for a stabilizing controller. Using 
a similar approach, the results of the paper can be extended to translation invariant interconnections of linear quantum 
stochastic systems on higher dimensional lattices. 
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